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ABSTRACT 

Forward shocks caused by the interaction between a relativistic blast wave and the circum- 
burst medium are thought to be responsible for the afterglow emission in Gamma-Ray Bursts 
(GRBs). We consider the hydrodynamics of a spherical relativistic blast wave expanding into 
the surrounding medium and we generalize the standard theory in order to account for several 
effects that are generally ignored. In particular, we consider the role of adiabatic and radiative 
losses on the hydrodynamical evolution of the shock, under the assumption that the cool- 
ing losses are fast. Our model can describe adiabatic, fully radiative and semi-radiative blast 
waves, and can describe the effects of a time-varying radiative efficiency. The equations we 
present are valid for arbitrary density profiles, and also for a circum-burst medium enriched 
with electron-positron pairs. The presence of pairs enhances the fraction of shock energy 
gained by the leptons, thus increasing the importance of radiative losses. Our model allows 
to study whether the high-energy (> 0.1 GeV) emission in GRBs may originate from after- 
glow radiation. In particular, it is suitable to test whether the fast decay of the high-energy 
light curve observed in several Fermi LAT GRBs can be ascribed to an initial radiative phase, 
followed by the standard adiabatic evolution. 

Key words: gamma-ray burst: general - X-rays: general 



1 INTRODUCTION 

Afterglow radiation in Gamma-Ray Bursts (GRBs) is attributed 
to external shocks produced by the interaction between the ultra- 
relativistic ejecta and the circum-burst medium (CBM). As the col- 
lision evolves, a relativistic forward shock propagates into the sur- 
rounding medium, sweeping up and accelerating the ambient mat- 
ter. Synchrotron emission from the shocked electrons powers the 
observed afterglow, which is usually detected in the X-ray, optical 
and sometimes radio bands. 

Recent observations by the Fermi/LAT (Large Area Tele- 
scope; lAtwood et al.l2009l) and AGILE/GRID (Gamma- Ray Imag- 
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120091 ; iBarbiellini et alj[2002h telescope 



ing Detector; 

have revealed that a few percent of GRBs, both long and short, 
are detected in the MeV-GeV energy range, confirming previous 
findings by EGRET (the Energetic Gamma Ray Experiment Tele- 
scope onboard the Compton Gamma-Ray Observatory). By in- 
specting the temporal and spectral behavior above 0. 1 GeV, some 
generic properties of the high-energy emission have been inferred 
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dGhisellini et alj20ld : lbmodei et al]2009l) : (i) the high-energy ra- 
diation typically starts at early times, during the prompt phase, with 
a delay up to a few seconds after the beginning of the sub-MeV 
prompt emission; (ii) it lasts longer than the prompt keV-MeV 
emission detected by the Fermi/GBM; (Hi) the high-energy flux de- 
cays as a power law in time, typically steeper than r 1 . 

The emission mechanism remains controversial. A promis- 
ing possibility, investigated by several authors, is that the high- 
energy emission has the same origin as the afterglow, at least when 
the MeV-GeV emission appears as a spectral component sepa- 
rate from the prompt, and it persis t s after the prompt has faded 
out dKumar & Barniol Durar] (20091 l20ld: iGhisellini et alj|20ld : 
[ Ghirlanda et al.ll2010l; ICorsi et al.ll20ld: |Gao et alj I2009L but see 
biran & Nakarll20ld : fMaxham et alj|201 1| for caveats). The bright- 
est bursts detected by LAT show a peak in their > 0. 1 GeV light 
curve at early times. If interpreted as the onset of the afterglow 
emission, the peak time can be used to estimate the initial bulk 
Lorentz factor r of the blast wave. The typical values obtained 
for the LAT bursts are quite large, around To~10 3 dGhisellini et al.l 
l20ld : lGhirlanda et al.l20l"d) . Beyond the peak, the > 0.1 GeV light 
curves decay with a temporal slope steeper than t , which can be 
easily explained if at these times radiative losses affect the dynam- 
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ical ev olution of the blast wave, as suggested by iGhisellini et all 

feoioh . 

Models where the energy lost to radiation is a sizable frac- 
tion of the blast wave energy are referred to as radiative. Radia- 
tive evolution sets in if the following two conditions are both satis- 
fied: (0 a considerable fraction of the energy dissipated in external 
shocks is given to leptons (e e - 1) and (ii) leptons are in the fast 
cooling regime, i.e. e rat i is close to unity. Here, e ra d is defined as 
the fraction of the energy gained by electrons which is promptly 
radiated. Estimates derived from modeling the X-ray and optical 
afterglow sugge st that e e m ay be roughly distributed around 0.1 
(Panaitescu & Kumar 2001), which has been often adopted as a 
typical value. Even if our knowledge of particle acceleration in rela- 
tivistic shocks is poor, it is reasonable to imagine that higher values 
of e e could be achieve d if the pre-shock med ium is enriched with 
electron-positron pairs. IGhisellini et"al, considered the pres- 

ence of pairs in the circum-burst medium to explain why the fire- 
ball may be radiative at early times, when the GeV emission takes 
place. A process that produces pairs in the ambient medium (invok- 
ing prompt photon s scattered by the circum-b urst electrons) has 
been c onsidered by Thompson & MadavJ d2000h and|Beloborodov 
j2002l) . This process is efficient up to some radius _/?i ad. beyond 
which the probability for an ambient electron to intercept a prompt 
photon becomes smaller than unity. 

Even when the fireball is highly radiative in its early evolution, 
a transition to the standard adiabatic regime should occur when: (/) 
the radiative process switches from fast to slow cooling (i.e. e md 
decreases from ~ 1 to e ra( /<K 1); or (ii) the physical conditions mo- 
tivating large e e cease to hold. For example, in the pair enrichment 
model described above, at radii larger than i? toad the pair loading 
gets inefficient, the pair-to-proton ratio drops and e e decreases, 
driving the transition from the radiative to the adiabatic regime. 

At present, nearly 30 GRBs have been detected by LAT and 
the relation between the high-energy emission and the afterglow 
radiation can now be investigated in more detail. To probe the pres- 
ence of an initially radiative epoch and to follow the transition to the 
adiabatic regime as the pairs-to-proton ratio in the CBM decreases, 
we need to generalize the standard model for the hydrodynamical 
evolut ion of GRB externa l shock s. 

iBlandford & McKed (1976) found a self-similar solution that 
describes the deceleration stage of adiabatic blast waves (i.e., in the 
case of negligible radiative losses), valid as long as the velocities 
are relativistic. Simpler analytical models based on the thin homo- 
geneous shell approximation (where the hydrodynamical proper- 
ties of the fluid behind the shock are taken to be uniform inside the 
shell) allow to describe the whole evolution of the shell Lorentz fac- 
tor, from the coasting to the non-relativistic epoch, but they suffer 
from some limitations. During the relativistic deceleration stage of 
adiabatic blast waves, they differ from the Blandford- McKee solu- 
tion by a numerical factor JPanaitescu & Kum ar 2000). Also, some 
important ingredients are often missing, such as the role of adia- 
batic expansion losses jDermer & Humi 2001) or a proper descrip- 
tion of the non-relativistic phase jHuang et al Jl 19991) . In all these 
approaches, the external medium density is usually considered to 
have a power-law radial profile and negligible pair loading. 

In this paper we develop a theoretical model able to describe: 

• the whole blast wave evolution, starting from the initial coast- 
ing phase (r = r = const) down to non-relativistic speeds; 

• adiabatic, fully radiative and semi-radiative regimes; 

• the transition from the radiative to the adiabatic regime, due to 
the change from fast to slow cooling or to the decrease in e e \ 



• adiabatic expansion losses; 

• a CBM characterized by an arbitrary density profile, possibly 
enriched with electron-positron pairs. 

Our mo del is b ased on the simple homogeneous shell approx- 
imation (see Piran 1999| for an exhaustive treatment of this ap- 
proach). However, we obtain a more realistic description (;') by 
recovering the Blandford & McKee (1976) self-similar solution 
during the relativistic deceleration phase of adiabatic blast waves; 
and (ii) by properly accounting for the adiabatic losses. By in- 
cluding the adiabatic losses in a self-consistent way, we can cap- 
ture the re-acceleration of the blast wave when the CBM den- 
sity drops faster than R 3 , and the internal energy accumulated at 
earlier times is converted back into bulk kine tic energy. Also, we 
correctly recover the Sedo v-Taylor solution dSedov|[l946l ; iTavioil 
ll95d : lvon Neumann|[l947h during the non-relativistic regime. 

Our model allows to investigate more general issues, not nec- 
essarily related to the GeV emission. In particular, in this work we 
use our model to study: 

• the importance of radiative losses for arbitrary values of e e ; 

• the impact of an initial radiative phase on the energy content 
of the blast wave, and its effect on the subsequent adiabatic epoch 
(and then on the X-ray and optical emission); 

• the relation between the initial bulk Lorentz factor r and the 
peak time of the afterglow light curve. 

In a follow-up paper (Nava et al., in preparation), we will test the 
external shock origin of the MeV-GeV emission, by simultane- 
ously fitting the GeV, X-ray and optical light curves and spectra 
of the GRBs detected by LAT. We will adopt the dynamical model 
developed in this work and we will consider both the synchrotron 
and the synchrotron self-Compton processes, in order to predict 
light curves and spectra at different frequencies. This study will 
probe whether an initial radiative phase is needed to explain the 
steep GeV light curves, and if the pair enrichment model described 
above can be a viable explanation for this early radiative epoch. 

The presence of such a radiative phase could also affect the 
time when the blast wave becomes non-relativistic (which should 
occur earlier if the fireball efficiently loses energy in its early evolu- 
tion). A detailed description of the transition to the non-relativistic 
phase, accounting for the effects of lateral spreading in the realis- 
tic case of a collimated outflow, will be the subject of a follow-up 
paper (Sironi et al., in preparation). 

In the following, primed quantities are measured in the frame 
comoving with the fireball, unprimed quantities are measured in 
the frame of the progenitor star. This will be true for all quantities 
except the particle momentum p and the particles random Lorentz 
factor y, that are comoving quantities, but that we leave unprimed 
for the ease of notation. 



2 AFTERGLOW EMISSION: THE STANDARD MODEL 

It is possible to identify different stages in the evolution of a spher- 
ical relativistic blast wave. After the acceleration phase, the ejecta 
freely expand with constant Lorentz factor To (coasting phase). 
When the internal energy of the shocked swept-up matter ap- 
proaches the explosion energy (i.e., T^mc 2 ~ E ), the blast wave 
enters a second stage, where both the shock wave and the shocked 
material decelerate, and the kinetic energy is gradually transformed 
into random internal energy and partially radiated. In the third stage 
of the evolution, the blast wave decelerates down to non-relativistic 
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velocities. Most of the detected afterglow radiation is emitted dur- 
ing the relativistic deceleration stage, via synchrotron cooling of 
the shocked electrons. 

In this section, we briefly revise the methods that are com- 
monly used to describe the evolution of external shocks in 
GRBs, and we emphasize the limitations of the different ap- 
pr oaches. In §2.1 we review t he self-similar solutions presented 
by Blandford & McKee ( 1976) for adiabatic blast waves and by 
ICohen et al.1 d 19981) for semi-radiative shells. Self-similar solutions 
can describe the complete radial profile of the hydrodynamical 
quantities inside the shell, but their validity is limited to the rel- 
ativistic deceleration phase. In contrast, in the so-called homoge- 
neous shell approximation, the whole blast wave evolution can be 
followed, from the coasting to the non-relativistic epoch, by con- 
sidering a series of inelastic collisions between the shell and the 
external medium and by solving the equations for the energy and 
momentum conservation of the system. However, the hydrodynam- 
ical properties of the fluid behind the shock are taken to be uniform 
inside the shell. In §2.2 and 2.3 we present two common approaches 
based on the homogeneous shell approximation. 



2.1 Blandford & McKee 1976 (BM76) 

iBlandford & McKed d 19761. hereafter BM76) presented a fluid dy- 
namical treatment of a relativistic spherical blast wave enclosed by 
a strong shock, under the assumption of adiabatic evolution. They 
considered the relativistic jump conditions across the shock and 
solved the continuity equations to derive the radial profile of the hy- 
drodynamic quantities inside the shell. They derived a self-similar 
solution describing the deceleration stage of an adiabatic relativis- 
tic shock wave propagating into a medium with a power-law den- 
sity profile p(R) = AqR~ s , where R is the distance of the shock 
from the center of the explosion. Their result represents the ultra- 
relativistic version of the Sedov-Taylor solution, which describes 
the non-relativistic stage of adiabatic shocks. 

By definition, in the adiabatic regime the blast wave energy is 
constant and equal to the initial value E . By integrating the energy 
density of the shocked fluid over the shell volume, BM76 found: 



12 -4s j 
17 -As 



-T l mc l = 



16nA c z 
17 -As' 



R T 



(1) 



where m is the ambient mass swept up at radius R, and T is the 
bulk Lorentz factor of the fluid just behind the shock. This relation 
prescribes the dependence of F on R, giving T oc R- (3 -- m . In the two 
common cases (homogeneous density, s = 0; and wind-like profile, 
s = 2), it gives T oc R~ 3/2 and T oc R~ l/2 , respectively. Even if this 
approach provides a reliable description of the deceleration epoch 
for adiabatic blast waves, both the coasting and the non-relativistic 
phases (and the transition between different stages of the evolution) 
cannot be properly described. 

BM76 also investigated the case of fully radiative blast waves. 
In this regime, they considered all the shocked matter to be con- 
centrated in a cold, thin and homogeneous shell. From energy and 
momentum conservation, they found that the evolution of the shell 
Lorentz factor T satisfies 



r- 1 



(m + M ) 2 (r + 1) 

M 2 (r - 1) 



(2) 



where M is the initial mass of the ejecta. Contrary to the BM76 
adiabatic solution in Eq. [T] this equation (based on the homoge- 
neous shell approximation) is also valid in the coasting epoch and 



for non-relativistic velocities. It is the ultra-relat ivistic generaliza- 
tion o f the momentum-conserving snowplow dOstriker & McKeg 
1988), characterized by a thin shell where all of the matter is 
concentrated, with an empty cold interior. As pointed out by 
ICohen et al.1 dl998h . the solution in Eq.[2]has been derived without 
accounting for the factor 4/3 that comes from the Lorentz transfor- 
mati on of the ene r gy den sity in a relativistically hot fluid. 

ICohen et al.l 0998) have looked for self-similar solutions to 
describe the relativistic deceleration stage of semi-radiative blast 
waves, i.e. the generalization of the BM76 solution (Eq.QJ, which 
only holds for adiabatic shells. They consider an adiabatic shock 
followed by a narrow radiative region and a self-similar adiabatic 
interior. In their approach, the shell cavity is hot and pressurized, 
as opposite to the empty and cold interior of the radiative solu- 
tion by BM76 in Eq.|2] In the fully radiative regime, their solution 
does not reduce to the momentum-conserving snowplow of BM76, 
but it rather represents the ultra-re lativistic generalization of the 
pressure-driven snowplow solution (Ostriker & McKee 1983). 



2.2 Piran 1999 (P99) 

IPiranI dl999l. hereafter P99) solved for the evolution of spheri- 
cal blast waves in the homogeneous shell approximation, with a 
generic radiative efficiency. The evolution of the blast wave Lorentz 
factor is derived by imposing energy and momentum conserva- 
tion in the collision between the blast wave and the external mat- 
ter. An analytic equation for T(R) (from the coasting to the non- 
relativistic phase) can be derived. This method suffers from some 
limitations. In the adiabatic regime, the predicted scaling law be- 
tween T and R during the relativistic deceleration stage is the same 
as deri ved by BM76, but the two solutions differ by a numerical 
factor ( Panai tescu & Kumad koOO). Moreover, adiabatic expansion 
losses are not considered. For non-relativistic velocities, Piran's 
solution gives /3 oc R \ As iHuanget all d 19991) pointed out, this 
scaling does not agree with the Sedov-Taylor solution (fi oc R~ 3/2 ), 
which is expected to describe the final non-relativistic stages of the 
fireball evolution. In the fully radiative case, P99 recovers the same 
solution derived by BM76 (Eq.|2j, which describes a momentum- 
conserving snowplow. 

2.3 Huang et al. 1999 (H99) 

In order to derive a relation between T and R valid also in the non- 
relativistic regime, iHuang etalJd 19991 hereafter H99) considered 
the homogeneous shell model of P99 but imposed a different equa- 
tion to describe the internal energy of the fireball. Their assumption 
states that the comoving internal energy of adiabatic blast waves 
should be E\ m = (T - 1) m c 2 Q 

This prescription for the internal energy mimics the role of 
adiabatic losses. The Rankine-Hugoniot jump conditions state that 
the average kinetic energy per unit mass is constant across a rel- 
ativistic shock (as measured in the post-shock frame), meaning 
that the random Lorentz factor of an element dm of shocked ma- 
terial equals the shell Lorentz factor T(m) at that time. At early 
times, when the shell Lorentz factor was larger, the material was 
shocked to hotter temperatures. It follows that, in the absence of 
adiabatic and radiative losses, the internal energy of the shell should 
be E' int = J [T(m) - I] dm c 2 . By imposing E' im = (T - 1) m c 2 , H99 

1 For fully radiative blast waves, E' jm = 0, and their solution coincides with 
the momentum-conserving snowplow of BM76 and P99. 
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implicitly mean that the effect of adiabatic losses is to keep the 
random Lorentz factor of all the swept-up matter equal to the cur- 
rent bulk Lorentz factor of the shell, regardless of the hotter tem- 
peratures of the material accreted at earlier times. Although this 
treatment of the adiabatic losses relies on an ad-hoc prescription, it 
allows to recover the correct scaling between ft and R predicted by 
the Sedov-Taylor solution, /? oc Rr 3 ! 2 . 

A more accurate treatment of the adiabatic losses, giving 
the correct scaling of f3(R) in the non-relativisti c regime without 
any ad -hoc assumption, has been considered by IPanaitescu et all 
( 1998) and lDermer&Humil ( l200ll) . but only for adiabatic shocks 
in power-law density profiles. In our model, that we describe be- 
low, we generalize their formalism to the case of semi-radiative 
blast waves propagating in a medium with an arbitrary density pro- 
file. 



3 OUR MODEL 

In this section, we propose a novel method to describe the evolu- 
tion of spherical relativistic blast waves, in the homogeneous shell 
approximation. Our model overcomes the problems and limitations 
present in other approaches, that we have summarized in ij2] In par- 
ticular, our method is able to describe: 

• the coasting, deceleration and non-relativistic phases, and the 
transition between these regimes. For adiabatic blast waves, the rel- 
ativistic deceleration phase reproduces the BM76 solution, and in 
the non-relativistic regime the shell velocity scales as ft oc R~ 3/2 , as 
predicted by the Sedov-Taylor solution. 

• adiabatic, semi-radiative and fully radiative blast waves. 

• the case of a time-varying radiative efficiency, either resulting 
from a change in the shock microphysics (i.e., in the fraction e e 
of shock energy transferred to the emitting leptons, or due to the 
transition from fast to slow cooling, which determines a decrease 
in the electron radiative efficiency e ra( /). 

• arbitrary density profiles (e.g., resulting from the interaction 
of the progenitor stellar wind with the surrounding interstellar 
medium), including the case of a circum-burst medium enriched 
with electron-positron pairs. 

• possible re-acceleration phases, caused by the conversion of 
part of the internal energy back into bulk kinetic energy. Re- 
acceleration is expected to occur if the blast wave encounters a 
sudden drop in the density of the external matter (i.e., the density 
profile becomes steeper than p oc R 3 ). 

We now derive the equations describing the evolution of the 
shell bulk Lorentz factor T as a function of the shock radius R. 

Total energy — The energy density in the progenitor frame can be 
obtained from the Lorentz transformation: e = (e' + P')T 2 - P', 
where the pressure P' is connected to the comoving energy density 
e' and the comoving mass density p' by the equation of state P' = 
(y - l)(e' - p'c 2 ). Here, y is the adiabatic index of the shocked 
plasma, which we parameterize as y = (4 + r _1 )/3 for the sake of 
simplicity, obtainin g the expecte d limits y = 4/3 for T » 1 and y =; 
5/3 for r — > 1 (see |Pe'enl2012l for a more accurate prescription). 
The total energy in the progenitor frame will be E = eV = eV'/T, 
where V is the shell volume in the progenitor frame, and V = TV in 
the comoving frame. The total energy of the shell in the progenitor 
frame is the sum of the bulk kinetic energy and the internal energy: 



where we have defined 
_ yT 2 - y + 1 

1 eff = -p , 



(4) 



to properly describe the Lorentz transformation of the internal en- 
ergy. Here, M + m = p'V is the sum of the ejecta mass M = 
E /T c 2 and of the swept-up mass m(R), and E'., = (e' —p'c 2 W ' is 
the comoving internal energy. As pointed out bv lPe'eJ d2012 ). the 
majority of current models use T instead of T e g, with an error up to 
a factor of 4/3 in the ultra-relativistic limit, when T e g — » (4/3)T. 

Formal equation for T(R) — The blast wave energy E m in Eq. [3] 
can change due to (f) the rest mass energy dm c 2 accreted from the 
CBM, and (if) the energy lost to radiation dE ra d = T e ffdE' md . The 
equation of energy conservation in the progenitor frame then reads 

d [r(M + m)c 2 + r eff E' jm ] = dm c 2 + T eff dE' rad . (5) 

The overall change in the comoving internal energy dE'- t results 
from the sum of three contributions: 



dE'„„ = dE' sh + dE' ad + dE' md 



(6) 



The first contribution, dE', = (T— l)dmc 2 , is the random kinetic 
energy produced at the shock as a result of the inelastic collision 
with an element dm of circum-burst materiaQ. The second term, 
dE' ad , is the energy lost due to adiabatic expansion. The adiabatic 
losses lead to a conversion of random energy back to bulk kinetic 
energy. This contribution is generally negligible, but it may become 
important when the CBM density decreases faster than p oc R~ 3 , as 



shown in 93.21 The third term, dE' md , accounts for radiative losses. 

From Eq. [5] we derive the equation for the evolution of the 
shell Lorentz factor: 



dT 

Jr 



( r e# + i)(r-i)c 2 f + r eff ^f 

(M + m)c 2 + E' ml ^ 



(7) 



T(M + m) c 2 + T eff E' int . 



(3) 



which reduces to Eq. 5 of iDermer &Humi| [2001) if we replace 
T e g with T. For a spherical shell, dm/dR = 4nR 2 p, where p(R) is 
the circum-burst mass density. The term T e ffdE' ad /dR, accounting 
for adiabatic losses, allows to describe the re-acceleration of the 
fireball when the density of the CBM is a fast decreasing function 
of R (see 93.2l and the Appendix). 

Specifying dE' ad and E' im — Eq. [7] is the formal equation giving 
r(K). To solve it, we need to spec ify dE'„ . and E'„ r . To this aim, we 
adopt the approach proposed by Dermer & Hum| d200ll) and we 
generalize their method to account for (f) arbitrary profiles of the 
circum-burst density p(R), and (fi) a time-varying overall radiative 
efficiency e(R). The efficiency e is defined as follows: a fraction e e 
of the energy dissipated by the shock dE' sh is gained by the leptons, 
which then radiate a fraction e m d of their internal energy. It follows 
that the energy lost to radiation is dE' md = -e mc ]e e dE' sh = —edE' sh , 
with e = e rac ie e . In other words, e is the overall fraction of the shock- 
dissipated energy that goes into radiation. 

After gaining a fraction e e of the shock energy, the mean ran- 
dom Lorentz factor of post-shock leptons becomes y acc e - 1 = (T — 
\)e e lp e . Here, fi e = pjp is the ratio between the mass density p e of 
shocked electrons and positrons (simply "electrons" from now on) 
and the total mass density of the shocked matter p. In the absence of 
electron-positron pairs, electrons and protons have the same num- 
ber density and the same radial profile, and p e is simply the ratio 

2 As pointed out by BM76, in the post-shock frame, the average kinetic 
energy per unit mass is constant across the shock, and equal to (r - 1) c 2 . 
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between the electron and proton mass: p e = m e /(m e +m p ) =i m e /m p . 
We assume that, right behind the shock, the freshly shocked elec- 
trons instantaneously radiate a fraction e ra d of their internal energy, 
so that their mean random Lorentz factor decreases down to 

J Tad* - 1 = (1 - £ ra d)(y a cc,e - 1) = (1 ~ e rad )(T - 1)— . (8) 

He 

Following this transient, electrons cool only due to adiabatic 
losses. The assumption of instantaneous radiative losses is verified 
in the fast cooling regime (e md ~ 1), which is required (but not 
sufficient) to have e ~ 1 (i.e., a fully radiative blast wave). In the 
opposite case e ra d <k 1, the evolution is nearly adiabatic (e <K 1), 
regardless of the value of e e , and the details of the radiative cooling 
processes are likely to be unimportant for the shell dynamics. The 
case with intermediate values of e ra d and e is harder to treat analyt- 
ically, since the electrons shocked at radius R may continue to emit 
copiously also at larger distances, affecting the blast wave dynam- 
ics. In the following, we implicitly assume either that the emitting 
electrons are in the fast cooling regime (e md ~ 1) or that the radia- 
tive losses are unimportant for the shell dynamics (i.e., e <k 1). 

Assuming that protons gain a fraction e p of the energy dissi- 
pated by the shock (with e p = 1 - e e - e B \ — e e , if the fraction e s 
converted into magnetic fields is negligible), their mean post-shock 
Lorentz factor will be 



r« lf , P -i = (r-D— 
h p 



(9) 



where p p = p p /p is the ratio between the mass density of shocked 
protons p p and the total shocked mass density p. In the standard 
case, when pairs are absent, p p & 1. Since the proton radiative 
losses are negligible, the shocked protons will lose their energy 
only due to adiabatic cooling. 

Comoving volume and adiabatic losses — To estimate the adia- 
batic losses, we assume that the shell comoving volume scales 
as V oc R 3 /T, corresponding to a shell thickness in the progeni- 
tor frame ~ R/T 2 . This scaling is correct for both relativistic and 
non-relativistic shocks, in the de celerating phas e (BM76). For re- 
accelerating relativistic shocks, IShapirol (l980) showed that the 
thickness of the region containing most of the blast wave energy 
is still ~ R/T 2 . For the sake of simplicity, we neglect changes in the 
comoving volume due to a time-varying adiabatic index or radia- 
tive efficiency. 

The radial change of the comoving momentum p of a shocked 
particle, as a result of expansion losses, will be governed by 



pdlnV _ p 1 3 d\nT \ 
~3 dr ~ ~3 \r ~ dr J 



(10) 



which is valid both in the relativistic and non-relativistic limits 
dDermer & Humil200ll) . The comoving momentum at radius R, for 
a particle injected with momentum p(r) when the shock radius was 
r, will be 



PaAR, r) = - 
K 



r(R) 



r(r) 



1/3 



p(r) 



(11) 



The comoving momentum p(r) is given by p = (y 2 - 1) where 
for electrons y = y ra d,e( r ) (Eq. [U, and for protons y = y acCtP (r) 
(Eq. [9). We remind that, by assuming that the shocked electrons 
promptly radiate at the shock, and then they evolve adiabatically, 
we are implicitly considering only the fast cooling regime e rai ~ 1 . 
Of course, our treatment is also valid in the quasi-adiabatic regime 
e« 1, when the radiative losses do not affect the shell dynamics. 



Internal energy and adiabatic losses — Considering the proton and 
lepton energy densities separately, the comoving internal energy at 
radius R will be 



E',JR) = 4ttc 



2 ^drr 2 {p p { 



r)-l]+p e (/)bw(«, r)- 1] (12) 



where y ac ] = {p 2 ai + 1) 1/2 . With the help of Eq.Qj] we can explicitly 
find E' jiu (R) and insert it in Eq.|7] 

The other term needed in Eq. [7] is dE' ad /dR. First, we derive 
(dyldR) al j for a single particle, using y = (p 2 + 1) 1/2 and Eq. 1101 
Then we integrate over the total number of particles, again consid- 
ering protons and leptons separately. Finally we have: 



dR 



\R 3 dR J Jo 



P P (r) 



ad,p 



3 dR 



+ Pe(r) 



(13) 



In Eqs.[T2]and[T3] we have assumed that only the swept-up matter 
is subject to adiabatic cooling, i.e., that the ejecta particles are cold. 

Relativistic phase — As long as the shocked particles remain rela- 
tivistic (i.e., y a d = p a[ ] » 1), the equations for the comoving inter- 
nal energy and for the adiabatic expansion losses assume simpler 
forms: 



EUR) = 4kc 2 drr 



f 

Jo 



r(r) 

r(R) 



1/3 



r(r) < —p p + (1 - e rad )—p, 

(Hp He 



1/3 



dR 



-EUR) 



T(r) 

i rfiogr 

3 dR 



T(r)p(r)(e 



(14) 
(15) 



In the second line of Eq. [14] we have used that e = e e e ra i. In the 
absence of significant magnetic field amplification, e p + e e ^ 1 so 
that e p + €g — e - 1 — e, and the radiative processes of the blast 
wave are entirely captured by the single efficiency parameter e. In 
the fast cooling regime e m d ~ 1 and e ^ e e . In this case the term 
€ p + €e — e reduces to e p , meaning that, regardless of the amount 
of energy gained by the electrons, in the fast cooling regime the 
adiabatic losses are dominated by the protons, since the electrons 
lose all their energy to radiation. 

Adiabatic regime — Evaluating these expressions for adiabatic 
blast waves in a power-law density profile p oc R \ we obtain 



EUR) = - — —Tmc 2 
"" y ' 9 -2s 



dE' ad (R) _ (9 - s)(3 - s) Tmc 2 
dR ~ 2(9 -2s) R 



(16) 



where we have used that T oc R {3 s)/2 as in the adiabatic BM76 
solution. 

Fully radiative regime — In the fully radiative regime e = 1, which 
implies E' im = and dE' ad = 0, Eq.|7]reduces to 



dT 

dm 



qy+ i)(T-i) 

M + m 



(17) 



which describes the evolution of a momentum-conserving (rather 
than pressure-driven) snowplow. If we replace T e g — > T, the solu- 
tion of this equation coincides with the result by BM76 (i.e., our 
Eq.0 and also with the equations derived by P99 and H99, since 
for fully radiative solutions the specific treatment of the adiabatic 
losses is unimportant. 
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Correction factor — Since our model is based on the homogeneous 
shell approximation, our adiabatic solution does not recover the 
correct normalization of the BM76 solution (see Eq.[T]l. In our treat- 
ment, the total energy of a relativistic decelerating adiabatic blast 
wave in a power-law density profile p(R) cc R~ s is 

E - T eff E' mt ^ ^TE' mt ^ ^^-T 2 mc 2 , (18) 

so that the BM76 normalization can be recovered if we multiply 
the density of external matter in Eqs. [7] [T2] and Q~3]by the factor 
(9-2^/(17-4^)0 To smoothly interpolate between the adiabatic 
regime (where we want to recover the BM76 solution) and the ra- 
diative regime, we suggest a correction factor: 

Cm^se + ^-fl-e). (19) 
17 - 4.? 

No analytic model exists that properly captures the transition be- 
tween an adiabatic relativistic blast wave and the momentum- 
conserving sn owplow, as e incre ases from zero to unity (instead, 
the model by ICohen et al.l 1998 approaches the pressure-driven 
snowplow in the limit e — » 1). For this reason, we decide to 
employ the simple interpolation in Eq. [19] in order to join the 
fully adiabatic BM76 solution with the fully radiative momentum- 
conserving snowplow. 

In summary, Eqs. [7] [12] and [13] complemented with the cor- 
rection in Eq. [19] (which should by applied to every occurrence of 
external density and external matter) completely determine the evo- 
lution of the shell Lorentz factor T as a function of the shock radius 
R. 

In our model it is easy to recover the solutions by P99 and 
H99. First of all, we should replace T e g — > T and Cbmi6,c 1. 
The approach described by P99 can be recovered for dE' ad = 0, so 
dE' jnt = (1 - e)(r - \)dmc 2 , while the equation given by H99 can 
be obtained for E' im = (1 - e)(T - i)mc 2 , which implies dE' ad = 
(1 - e)dTmc 1 , at odds with respect to our Eq. 1 131 




R/R 

' dec 

Figure 1. Evolution of the shell momentum r/J as a function of the shock 
radius R, normalized to the deceleration radius Rj ec , in the case of homoge- 
neous ISM density (s = 0, upper panel) and wind-like profile (.v = 2, bottom 
panel). Solid curves show the prediction of our model for the fully adiabatic 
and fully radiative solutions. Dotted (dashed) curves show instead the P99 
(H99) solution. The straight (yellow) line is the prediction by BM76, which 
is valid after the deceleration started and till the velocities are relativistic. 
The grey area indicates the non-relativistic regime. 



3.1 Testing the model 

Before applying our model to the physics of GRB afterglows, it is 
healthy to compare it with previous studies and to understand the 
differences. 

Momentum of the blast wave: adiabatic and fully radiative cases 
Fig. H shows the evolution of the shell momentum F/J for a CBM 
density profile p = AqR~'. We compare our model (solid curves) 
with the predictions by P99 (dotted) and H99 (dashed), for adia- 
batic and fully radiative blast waves, for homogeneous (s = 0, top 
panel) or wind (s = 2, bottom panel) density profiles. The shock ra- 
dius is normalized to the deceleration radius, which is defined such 
that the swept-up mass is mj ec = m(Rd ec ) = Mq/Tq, so 

l/(3-s) 



R.i. 



(3 - s)M 



47iA n r n 



(20) 



In the relativistic adiabatic regime, Fig. [T] shows that our solu- 
tion reproduces the adiabatic relativistic solution of BM76 (yellow 
straight lines) for both the uniform and the wind density profile. 
P99 and H99 underestimate the shell momentum with respect to 
the BM76 solution. Moreover, the approach by P99 does not cor- 
rectly reproduce the Sedov-Taylor solution (with fj oc R~ 3/2 ), since 
it gives f3 oc R 3 (as anticipated in 5j2). For fully radiative blast 

3 This correction is only appropriate for a power-law density profile. 



waves, the solutions by P99 and H99 coincide and differ from our 
solution only because they use T instead of T e ff in Eq .1171 

As long as we are in the relativistic regime, Fig. Q] shows that 
all models agree on the slope of Tp as a function of R, both for the 
adiabatic and for the fully radiative case. 

Momentum of the blast wave: semi-radiative case — Fig. [2] shows 
the slopes a = -dlnT/dlnR as a function of the overall radiative 
efficiency e. All three methods give a = (3 - s)/2 for e = and 
a = 3 - s for e = 1, but they differ for intermediate values of 
6 (i.e., they predict different slopes in the case of semi-radiative 
blast waves). In this case the value of a in the model by P99 can 
be derived analytically, and it gives a = (3 - s)/(2 - e). In our 
model, an analytic expression for a can only be obtained for e <k 1, 
and it yields a = (9 - 3s + 3e)/(6 - 4e). The same is true for 
the approach by H99, which results in a = (3 - s)/(2 - 2e) for 
f « 1. From Fig.[2] it is clear that the model by P99, which neglects 
the adiabatic expansion losses, always gives the flattest slopes (i.e., 
the shell evolution remains close to the adiabatic case, everything 
else being fixed). Comparing our model with H99 we see that our 
model yields flatter slopes for a homogeneous medium (top panel), 
whereas the opposite holds for a wind profile (bottom panel). 

Partition of energy — Fig. [3] shows the partition of energy among 
different components (red for kinetic energy, green for internal en- 
ergy, purple for radiated energy), during the evolution of the blast 
wave. The results of our model are plotted with solid and long- 
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Figure 2. Comparison of the slope a (such that T oc R~ a ) in our model 
(solid lines) with the predictions by P99 (dotted) and H99 (dashed), for 
different values of the radiative efficiency e. The slope is measured during 
the relativistic deceleration phase of evolution. The top panel refers to the 
case of a homogeneous interstellar medium, the bottom panel to a wind 
density profile. We remark that the slopes in our model do not depend on 
the correction factor Cbm16.c m Eg. 1191 

dashed lines, whereas the predictions by P99 and H99 are shown 
with dotted and dashed lines, respectively. In the coasting phase 
(R < Rd ec ), the blast wave energy is still dominated by the bulk 
kinetic energy of the ejecta (red curves). At R ~ R c i ec , bulk and 
internal energies are comparable by definition. For R > Rd ec , the 
ejecta decelerate, and the contribution of the bulk motion to the 
overall energy budget decreases. 

For adiabatic blast waves (top panel), the decrease in kinetic 
energy (red lines) is balanced by an increase in internal energy 
(green lines). After the shell has accreted a mass m ~ M , the blast 
wave kinetic energy becomes dominated by the swept-up material, 
rather than by the initial ejecta mass. This leads E kj „ to increase as 
E kin ~ Tmc 2 oc m l/2 . In the non-relativistic regime (shaded area), 
the internal energy E im = T(r - l)m c 2 ~ (r - l)m c 2 and the bulk 
kinetic energy E kj „ = (T - 1)(M + m)c 2 ~ (T - \)mc 2 asymptote 
to the same value in the model by H99 (dashed lines). In our model 
the internal energy dominates over the kinetic energy for a uniform 
medium (solid lines), whereas the opposite holds for a wind-like 
density profile (long-dashed lines). Finally, we remark that in the 
model by P99 (dotted lines), that neglects the adiabatic expansion 
losses, most of the energy at late times is still in the form of internal 
energy, with a negligible contribution by the bulk kinetic energy. 

For radiative blast waves (bottom panel), the kinetic energy 
(red line) is lost to radiation (purple line). The blast wave is cold 
(i.e., the contribution by the internal energy is negligible), since the 
random energy generated at the shock is promptly radiated. 

3.2 Complex density profiles and re-acceleration 

The formalism described in this paper is valid for any spheri- 
cally symmetric profile of the external density. If the blast wave 
encounters a change in the CBM profile during the deceleration 
phase, with density decreasing faster than p oc R 3 beyond some 
radius, the conversion of the accumulated internal energy back 
into bulk kinetic energy can lead to a re-acceleration of the blast 



TOTAL ENERGY £ = 




1 10 1 2 

(R/R d J (3 " s)/3 

Figure 3. Evolution of the different contributions to the blast wave energy, 
normalized to the initial energy Eg. Solid (long dashed) lines refer to our 
model in the case of s = (.v = 2). Dotted and dashed curves refer to 
P99 and H99, respectively. The top panel shows the adiabatic case e = 0, 
the bottom panel presents the fully radiative case e = 1 . The contribution 
of the bulk kinetic energy £t,„ = (r - 1)(Mq + m)c 2 is shown in red, the 
internal energy E{ nl = T e gE' jnl m green, the energy lost to radiation E m j = 
J T e jj dE' md in purple. The sum of the different contributions is the cyan 
curve, demonstrating the conservation of energy. The grey area indicates 
the non-relativistic regime. 



wave. As shown by IShapird dl980h . the Lorentz factor of re- 
accelerating blast waves should scale as T oc /j( J - 3 >/ 2 . of course, 
the re-acceleration happens only for adiabatic blast waves; if the 
shell is fully radiative, all the internal energy is promptly radiated 
away, and it cannot produce any re-acceleration. 

We consider the case of an adiabatic blast wave propagating in 
a medium where the external density is constant up to R b = 15/?^ 
and then it decreases as p oc R~ 1/2 . Fig.|4]shows the evolution of the 
shell momentum. Our model (red solid line) shows that the evolu- 
tion of the shell momentum departs from the result expected for a 
homogeneous medium (red dashed line) at radii R > Rb, where the 
shell starts to re-accelerate. Its momentum increases as T oc R [/4 
(red dotted line), as expected for a density profile with s = 7/2. 
At large radii, the shell Lorentz factor asymptotes to a value which 
is lower than T (we find T ~ 10 2 in the example in Fig.|4j. This 
is just a consequence of energy conservation. Only the internal en- 
ergy stored in the shell at the end of the deceleration phase (R = Rb) 
can be used to re-accelerate the blast wave. If m b = m(R b ) is the 
mass accreted at R b , the asymptotic value of the shell Lorentz factor 
can be obtained directly from the conservation of the bulk kinetic 
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Figure 4. Evolution of the momentum T/3 for an adiabatic blast wave prop- 
agating in a structured density profile, which passes from a uniform medium 
to ap cc R~ 7 ' 2 profile at Rt = 15Rj ec . The solution of our model is shown 
with the red solid line, whereas the predictions by P99 and H99 are plotted 
with a dotted and a dashed line, respectively. It is apparent that the mod- 
els by P99 and H99 cannot capture the re-acceleration of the blast wave 
occurring at R > 15 Rj ec . We also show the solution corresponding to uni- 
form d ensity every where (red dashed line), the scaling T oc R 1 ' 4 predicted 
by Shapiro 1 1980) for the re-acceleration phase (red dotted line), and the 
expected saturation at late times (red dot-dashed line). 



energy: 



r-i = (r -i)- 



M + ■ 



(21) 



-m, 



where the factor (9 - 2i)/(17 - 4s) comes from the correction 
we adopt to recover the normalization of BM76 for adiabatic de- 
celerating blast waves (see Eq. |19t. The capability to capture the 
blast wave re-acceleration is peculiar to our model and it requires a 
proper description of the adiabatic expansion losses. In contrast, the 
methods by P99 or H99 cannot describe the re-acceleration phase 
(dashed and dotted black curves in Fig. [4] respectively). 

We find that the evolution of T with radius in the re- 
accele rating phase departs from the scaling predicted by Shapiro 
dl980h for large values of s, and it approaches T oc R 3/4 in the limit 
s » 1 , which is the expected behavior for a fireball expanding in 
the vacuum, as we demonstrate in Appendix [A] 

In summary, our model can be used to follow the dynam- 
ics of blast waves that propagate in a structured CBM, with den- 
sity enhancements (or "clumps") and troughs, as it is often con- 
sidered to explain flares and rebrightening s in the optical and X- 
ray light curves of many GRB aftergl ows ( Nak ar & Granoll2007l ; 
iMesler et alj|201 j : lHascoet et alj|2012h . However, we remark that 
our results, being derived in the homogeneous shell approxima- 
tion, cannot describe properly the role played by reverse shocks 
and rarefaction waves created when the shell interacts with a sud- 
den density jum p or drop. For this, on e needs specific analytical 
treat ments (e.g., Nakar &Granot| [2007) or numerical simulations 
(e.g. jMimica & GianniosfeOl 1 ). Also, in the present work we only 
consider a spherical shell, and we defer the discussion of jetted out- 
flows (and their associated lateral expansion) to a follow-up paper 
(Sironi et al., in preparation). 



4 RESULTS 

4.1 Bolometric light curves 

To predict the afterglow light curves we need to estimate how the 
observed luminosity evolves as a function of the observer's time. 
We only focus on the bolometric light curve, so we do not need 
to specify the emission mechanism. We proceed in two steps: first, 
we estimate the luminosity as a function of the shock radius R, and 
then we introduce the relation between R and the observer's time t. 

The luminosity observed from a parcel of shocked fluid mov- 
ing at an angle 9 with respect to the line of sight will be L g (9) = 
6 4 L', where 6 = [T(l -/Jcosf)] -1 is the Doppler factor of the 
shocked fluid, and L' is the comoving luminosity. Assuming V to 
be isotropic in the comoving frame, the observed (isotropic equiv- 
alent) luminosity L is calculated by averaging S A over the angle 9 



alent) luminosity L is c 
( Ghirl anda et alj2012t) : 



L' 

T 



f 

Jo 



^ sin 9 d9 = T 2 ^4^-L' 



(22) 



In our formalism the comoving luminosity is L' = dE' md /df, where 
dE' md = e(T — \)dmc 2 is the energy lost to radiation in the fluid 
comoving frame, and for spherical explosions dm = 4nR 2 pdR. By 
employing this expression for dE' mi , we are implicitly saying that 
the emission at any given time only comes from the mass element 
dm just passed through the shock. In other words, we are assum- 
ing that the electrons radiate right behind the shock, and then they 
evolve adiabatically. It follows that the expression we find below 
for the observed luminosity is valid only in the fast-cooling regime. 

We then need to compute the distance dR traveled by the shock 
during a comoving time interval df, which turns out to be 



dR : 



A/ 



-cdf , 



(23) 



r(i -pp A ) 

where fi sh is the shock velocity in the progenitor framed This ex- 
pression can be simplified by employing the convenient parame- 
terization /?,/, = 4r 2 /?/(4r 2 - 1), which is valid both in the ultra- 
relativistic phase (giving T 2 sh =; 2T 2 ) and in the non-relativistic 
limit (fi sh m 4/3/3). In this relation for /?,/,, the speed /? should 
be interpreted as the fluid velocity right behind the shock, before 
the radiative cooling losses have appreciably chan ged the dynam- 
ics of the post-shock flow (see lCohen et al.ll l998b- It follows that 
(1 - PPsh) 1 = r 2 (/3 2 + 3)/3, and the comoving luminosity is 

L ' = e(r -i )r ^±l^ c 3 , (24) 

3 dR 

so that the observed bolometric luminosity as a function of the 
shock radius will be 



L = e(r- 1)T 3 



/3 2 + 3 



dm , 
P «-R C 



(25) 



The linear dependence on the shock velocity is expected, since 
in our model the luminosity is proportional to the rate at which 
the shock picks up matter from the circumburst medium. In the 
non-relativistic limit, the observed luminosity reduces to L = 
6(l/2y3 2 )y6,;, dm/dRc 3 , as expected from a shock with velocity f} sh 



4 One could naively expect that dR = fi s i,Tcdf . The difference with respect 
to the expression we employ comes from the fact that, in our approach, a 
given fluid element radiates only right after passing through the shock, so 
two subsequent emitting events are spatially coincident in the shock frame, 
rather than in the fluid comoving frame (as implicit in dR = /3 s i,Tcdt'). 
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Figure 5. Blast wave energy content (normalized to its initial value) as a 
function of the observer's time for a blast wave with constant efficiency e. 
Different values of e (from 0.01 to 1) are considered, both for .s = (solid 
lines) and s = 2 (dashed lines). For comparison, the solution proposed by 
Islrilll997l) for e = 0.1 and s = is also shown (dotted black line). 




Figure 6. Blast wave momentum T/3 as a function of R for a shell with ra- 
diative efficiency e that evolves in radius as shown in the inset. Two differ- 
ent cases for e(S) are considered (purple and green lines). The dashed grey 
curves refer to the standard adiabatic and radiative cases with e constant 
and equal to 10~ 2 and 1, respectively. The dotted and dot-dashed purple 
curves show the adiabatic solutions with e = const = 10~ 2 having the same 
late-time energy as the purple solid curve, and either the same To (purple 
dotted line) or the same Mq (purple dot-dashed line). 



that converts a fraction e of the post-shock internal energy per unit 
mass (1/2) f3 2 c 2 into radiation. 

Finally, to obtain the luminosity as a function of the observer's 
time, we need to relate the shock radius R to the observer's time t. 
For the sake of simplicity, we do not consider the formali sm of 
equitemp oral surfaces dis cussed by iBianco & Ruffini j2005l) . but 
following lwaxmarj i 19971) . we use the relation 



t(R) = t R + t e - 



Jo / 



dr + ■ 



R 



(26) 



p sh c " F-(l+/3)c 

which consists of the sum of a radial delay t R , i.e., the difference 
between light travel time to radius R and shock expansion time to 
the same radius, and an angular delay t s , for photons emitted from 
the same radius R but at different angles with respect to the line of 
sight. The angular term t s dominates in the ultra-relativistic limit, 
whereas the radial term t R prevails in the non-relativistic phase. 

We point out that the expression for the observed luminos- 
ity in Eq. [25] and the relation t(R) in Eq. [26] are valid both in the 
ultra-relativistic and non-relativistic phases. From their combina- 
tion, one could estimate the asymptotic scaling of the observed lu- 
minosity with time, as the blast wave decelerates. For adiabatic 
blast waves, we find L oc £(t)E r l both in the relativistic and 
non-relativistic regimes (provided e(f) <K 1 at all times). Note 
that in both regimes the bolometric luminosity does not depend 
on the value of the CBM density. For radiative blast waves we 
find L oc [M*- 2, r <10 - M /A ] 1/<7 ~ 2l) in the relativistic phase and 
L oc [M^ , r <7_2s) /'4o] 1/<4 in the non-relativistic limit (assuming 
e(t) ~ 1 at all times). 



4.2 From the radiative to the adiabatic regime 

One of the crucial parameters that could help to unveil the nature 
of the progenitor of GRBs is the total energy output of the central 
engine E tot . This quantity can be estimated from the sum of the 
energy emitted during the prompt phase (E y ) and the energy of the 
blast wave that powers the afterglow emission (E ). Moreover, from 
these two quantities one can derive the efficiency of energy dissipa- 
tion in the prompt phase rj = E y IE tou which could help to constrain 
the uncertain mechanism that produces the prompt y-rays. 

Estimates of Eq require observations of the afterglow over 



a broad range in time and frequency (e.g., Panaite scu & Kurriail 
2002; Yost et alj|2003l) . An alternative method is based on obser- 
vations at frequencies above the synchrotron cooling frequency v c , 
since they give a robust proxy for the tot al blast wave energy E 
dKumarl2000l : lFreedman & Waxmanl200ll) . For typical parameters, 
the X-ray data at late times lie above v c , so the X-ray flux can be 
used to constrain the blast wave energy E and the efficiency of 
the prompt emission rj = E 7 /(E y + Eq). Several authors have ap- 
plied this method to different samples of both long and short GRBs 
dLloyd-Ronning & Zhanel 120041 : IZhang et al J 120071 : iRacusin et all 
201 1 ; Berger 2007), finding similar results: rj reaches values of 10% 
or higher, challenging the simplest version of the internal shock 
model, where only a few per cent of th e jet energy can be dissipated 
( lKumar[[l999l : |Panaitescu et al.ll 19991) . Swift observations make the 
situation even worse, r aising some of the inferred GRB prompt ef- 
ficiencies up to >90% jZhang et alj2007l) . 

The problem can be alleviated if the afterglow energetics is un- 
derestimated. The adiabatic model cannot be strictly correct, since 
radiation takes away some energy from the blast wave. The X-ray 
flux is a proxy for the energy content of the blast wave at late times, 
but this is smaller than the initial energy Eq, if the evolution is 
not perfectly adiabatic. By neglecting the radiative losses, one may 
underestimate E and then overest imate r). Ra diative corrections 
based o n the equation proposed b y Sari { 19971) have been consid- 
ered bv lLlovd-Ronning & Zhand d2004l) . Even in this case, the in- 
ferred efficiencies rj are relatively high (between 4 0% and 100%). 
The validity of the equation derived bv lSaril dl997h has been ques- 
tioned, since it employs the relation R = 16r 2 cf between shock ra- 
dius anrlob^rver^s time, instead of the correct expression in Eq.|26l 
(seelWaxman 1997). So, the importance of radiative corrections in 
the afterglow energetics at late times is still unclear. In Fig. [5] we 
use our model to estimate, for constant values of e ranging from 
10~ 2 to 1, the energy content of the fireball E(t) as a function of 
the observer's time, normalized to t = R dec /(TgC). As an example, 
for a fireball with E = 10 53 ergs and r = 10 3 propagating in a 
medium with uniform density «o = 1 cirT 3 , the deceleration time 
is t ~ Is. The plot in Fig. [5] shows that, for a constant radiative 
efficiency e = 0.1, the energy content of the blast wave decreases 
by a factor of two in the first 10 4 seconds. For higher values of e, 
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Figure 7. Solid line: light curve of a blast wave with time-dependent e. 
The evolution of e(t) is shown in the inset. The grey dashed line shows the 
light curve of a blast wave with the same values of Eq and To, but e = const 
and equal to the value reached by the time-dependent e(t) at late times (i.e., 
e = 0.1). Dot-dashed and dotted lines show the light curves of two blast 
waves with constant e = 0.1, that have an initial energy equal to the late 
time energy content of the blast wave with time-dependent e. In particular, 
the dotted curve refers to the case of identical initial Lorentz factor To, while 
the dot-dashed line refers to the case of identical initial mass Mq . 



the fraction of fireball energy lost to radiation will be even larger . 
For comparison, we also plot the solution proposed bv lSaril fl997h 
for e = 0.1, which tends to overestimate the radiative losses. 

The results presented in Fig. [5] rely on the assumption that 
the efficiency e remains constant in time. The results significantly 
change if GRB afterglows experience an early fully-radiative phase 
followed by the standard adiabatic evolution at later times. To ad- 
dress this case, in Fig. [6] we study the evolution of the blast wave 
momentum TfS when e decreases from e = 1 to e = 1CT 2 (green 
and purple solid lines), and we compare it to the cases with con- 
stant e = 1 and constant e = 1CT 2 (dashed lines). We consider two 
possible choices for the time evolution of e, as shown in the inset 
of Fig. [6] (0 a step function where e = 1 up to a few R iec and 
e = 1(T 2 at larger radii (purple line); or (/() e = 1 for R < Rj ec and 
£ oc Rr - 1 for R > Rj ec (green line). In the latter case (green solid 
line), the blast wave momentum TyS initially follows the radiative 
solution (dashed line with £ = const =1), and then smoothly ap- 
proaches the scaling R~ i/2 expected for adiabatic blast waves. At 
large radii, the difference in normalization with respect to the so- 
lution having e = const = 10~ 2 (dashed line) comes from the fact 
that the fireball with evolving e has lost some of its energy during 
the early radiative phase. 

Consider now the case of a step function for the evolution of 
£(R) (purple solid line in Fig.[6]l. The blast wave momentum ini- 
tially follows the radiative solution, but at radii larger than 4i?^ c we 
have £(R) = 10~ 2 and the energy content remains nearly constant 
E(R > 4Rd ec ) ~ const = Eqj. At large radii, the blast wave evolves 
as a fireball with £ = const = 1(T 2 and initial energy Eq - Eqj. 
This is demonstrated in Fig.|6]for two different cases: (0 a fireball 
with Eq = Eqj and the same ejecta mass Mq as the shell with evolv- 
ing £ (purple dot-dashed curve); and (ii) the case of a fireball with 
E = E f and the same initial Lorentz factor r as the shell with 
evolving e (purple dotted curve). 

Let us now compute the bolometric light curve in the case of a 
time-evolving e. We are interested in reproducing the observations 
at GeV energies in the scenario of an early radiative phase of the 
GRB external shock. At the same time, we want to be consistent 
with the standard adiabatic evolution at later times. The light curve 



in the initial radiative regime (lasting a few tens of seconds) will be 
characterized by a steep temporal decay s imilar to that observed in 
the GeV emission dGhisellini et alj|201(j|) . which will be followed 
by the standard adiabatic scaling Leer 1 seen in the X-ray and op- 
tical bands. In this simple scenario we do not attempt to explain the 
plateaus and flares often present in the X-ray and/or optical data. 
Although bolometric light curves cannot be directly compared to 
the observations in a given frequency band, they are still useful to 
estimate how much energy can be radiated and lost by the blast 
wave before the adiabatic regime sets in. Consider the evolution 
of £ shown in the inset of Fig. [7j At early times, the blast wave is 
fully radiative (i.e., e(f) = 1). After a few seconds, the radiative 
efficiency decreases, reaching at t ~ 10 2 s the final value £ = 0. 1, 
which is the typical value inferred from late afterglow data. The re- 
sulting light curve is shown with the solid curve in Fig. [7] We also 
present, for comparison, the light curve of a fireball with the same 
initial values of Eq and To, but with e = const = 0.1 (dashed grey 
line). After 10 2 s, the two light curves have the same slope L oc r [ , 
and the difference in their normalizations gives the ratio between 
the initial energy Eq and the late-time energy content of the fireball 
with time- varying e. In the case shown in Fig.|7J this ratio is ~6, im- 
plying that ~ 83% of the initial energy Eq has been radiated in the 
first 10 2 seconds. By using observations after 10 2 s, as it is usually 
the case for optical and X-ray data, one would underestimate the 
blast wave energy Eq, and so overestimate the prompt efficiency n. 
Broadband fits of GeV, X-ray and optical data are needed in order 
to understand if an early radiative phase is required to fit the GeV 
light curves and, if this is the case, to estimate how much energy 
is radiated before the afterglow brightens in the X-ray and optical 
bands (Nava et al., in preparation). 



4.3 Estimate of the initial bulk Lorentz factor r 

During the initial coasting phase, when the swept-up matter is neg- 
ligible and the shell Lorentz factor is constant, the bolometric lu- 
minosity scales with time as f 2 ~ ! (i.e., it is rising for s < 2). When 
the deceleration starts (approximately at Rj ec ), the Lorentz factor 
decreases, producing a feature in the light curve: a peak (if s < 2) 
or a change in the power-law decay slope (if s > 2). The transition 
between the two regime s occurs approximately at t = Rdecl(T 2 c ) 
(Rees & Meszaros 1992), where R dec has been defined in Eq. |50] 
Since fo depends on the initial Lorentz factor (and it is nearly insen- 
sitive to all other parameters), the measurement of this time can be 
used to infer r dRees & Meszaro s 1992; Sari & Piran 1999). Since 
the peak time t peak of the afterglow light curve is a good proxy for 
Iq, we focus in this section on the relation between the peak time 
and r . We consider both adiabatic and semi-radiative regimes. 

Adiabatic case — Both before and after the peak time t p< , a u, the lu- 
minosity scales in time as a power law. From Eq. [25] by setting 
T = r we derive L(t < t peak ), while using T(R) given by the BM76 
solution we derive L(t > t pea k) (in both epochs we also use the 
relativ istic approximation /? ~ f} sh ~ 1). Following Gh irlanda et al.l 
l l2012l) . we extrapolate the two power-law segments of L(t) and find 
their intersection time, t A . The initial bulk Lorentz factor Tq is re- 
lated to the intersection time f A by the equation: 



r„(e = 0, s) 



(17-4.s)(9-2i)3 2 



110-2.! 



(4 - s)tt c 



5-j 



Agt 



3-s 



(27) 



Even if t A does not coincide with t pmk , for adiabatic models and 
s = it is possible to find numerically that t A = 0.9 t peak . Based on 
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our model, the relation between t peak and T for an adiabatic blast 
wave expanding in a homogeneous medium is then: 



r (e = 0, s = 0) = 190 



V,53 



peak.2 



Va.2 



1/8 



(28) 



where £ y-53 = £ r /10 53 ergs, t/ .2 = >7/0.2, and r pMiU = t peak / 100 s. 
We have replaced £0 by considering the relation E ^ £ y /j/, where 
£ r is the energy emitted in the prompt phase and rj is the radiative 
efficiency of internal shocks. We now compare our formula with 
those proposed in the literature for the adiabatic regime and s = 0. 
The equations m ost frequently used to estimate To fr om t mak are 
those derived bvlMolinari et al. 2007 (M07 hereafter) llLiang et al.l 
l20ld:lGruber et alj201 ll ; lRvkoff et al.l2009l ; lEvans et alj201 lh and 
bv lSari & Piranll 19991 (SP99). SP99 suggest the formula: 
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(29) 



while for M07 this expression represents Tj ec (the Lorentz factor 
at the deceleration radius Rj ec ), and to derive T they arbitrarily 
assume r = 2T dec . Note that our formula differs from the existing 
ones only by a numerical factor, since the dependences of To on Eq, 
A , and t peak are exactly the same. We estimate both numerically 
and analytically the discrepancies between our formula and those 
proposed by M07 and SP99. We find that our estimate of To is a 
factor of 1.64 systematically lower than the estimate by M07 and a 
factor of 1.20 larger than the estimate by SP99. These numbers are 
quite insensitive to the choice of Eq and Aq. 

We conclude that the equation proposed by M07 leads to a 
significant overestimate of T (by a factor of 1.64), while the equa- 
tion presented by SP99 gives a lower estimate of To, in better 
agreement with our result. The main source of disagreement be- 
tween our model and the approach by M07 is in their definition 
of the deceleration radius. In their model, the peak time t pea k of 
the light curve is assumed to correspond to the deceleration time, 
when the fireball reaches the deceleration radius. While SP99 de- 
fine the deceleration radius as the radius where the swept-up matter 
is m(R dec ) = M /T = E /(TgC 2 ) (i.e., the same definition adopted in 
this work), M0 7 use the relation m(R decM oi) = Eol{T dec c 2 ) (see also 
iLii et al .120121) . Fig.[8]shows that the luminosity peaks close to Rd ec , 
which is then a good proxy for the onset of the deceleration epoch. 
The deceleration radius Rjecmoi defined by M07, instead, overesti- 
mates the radius at which the luminosity reaches its maximum by a 
factor of ~ 2. This difference becomes even larger when consider- 
ing the light curve derived with the P99 model (dashed blue line). 

Semi-radiative case — The same approach can be applied to 
generic values of e, but in this case it is not possible to find an an- 
alytic formula that describes the Lorentz factor after the peak time 
(and then an analytic equation for L(t > t peak )). We then define: 



r„(e, s) = K(e, s) 



Ep 



K\e, s) 



(30) 



peak J 



and we list in Tab.[TJthe values of K and K' numerically estimated 
for different choices of e and s. 



5 CONCLUSIONS 

We consider the interaction of a spherical relativistic blast wave 
with the surrounding medium and we describe its dynamical evo- 
lution by taking into account radiative and adiabatic losses. We 



Table 1. Normalization factors for the estimate of To as a function of t A 
(K) and t peak (K'), calculated for different values of e and s (see Eq. |30t . 
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Figure 8. Light curve as a function of R/Rdec (where Rdec is defined as in 
Eq.|20}. The red solid curve is obtained by using our model, while the blue 
dashed curve is based on the P99 model. The vertical dashed and dotted 
lines mark the peak of the light curves for the two models. The dot-dashed 
line marks, instead, the deceleration radius as defined by M07. 



derive the equation that governs the evolution of the blast wave 
Lorentz factor as a function of the shock radius (Eq.[7J. This equa- 
tion is general and independent of the composition of the external 
medium and of the nature and efficiency of the emission processes. 
To present a simple solution of this equation, however, we restrict 
ourselves to the case where radiative losses are important only for 
the freshly shocked matter (which is strictly true only in the fast 
cooling regime), or to the case where radiative losses are not im- 
portant at all (quasi-adiabatic regime). The equations we derive 
are valid for any generic (spherically symmetric) density profile. 
We also give the correction factor that should be employed, in our 
homogeneous thin shell approximation, to match the Blandford & 
McKee (1976) solution in the relativistic deceleration stage. 

• Radiative losses — Our model is valid for any value of the 
overall radiative efficiency e = e e e rall , where e e , as usual, is the 
fraction of shock-dissipated energy gained by the electrons, while 
e ma - is the fraction of this energy which is lost to radiation. The 
value of e can evolve with time, thus varying the impact of the 
radiative energy losses on the shell dynamics. The decrease in the 
blast wave energy content due to the radiative losses is taken into 
account self-consistently in our approach. 

We study models where the value of e is kept constant in time and 
we derive, for each value of 6, the evolution with time of the energy 
content of the blast wave (Fig. [5). For e>0.1, the radiative losses 
can decrease the energy content of the fireball within t ~ 10 4 s by 
more than a factor of two. In the extreme limit e— 1, the energy 
content decreases by an order of magnitude in the first ~10 2 s. As 
shown in Fig. [2] radiative losses change the slope a of the scaling 
in radius of the blast wave Lorentz factor (such that T(R) oc R~ a ), 
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from the classical a = (3 - s)/2 of the adiabatic BM76 solution to 
the steeper a = (3 - s) of the fully radiative solution. 

We have also considered the realistic case where e varies in time 
from e =s 1 (fully-radiative) down tof « 1 (quasi-adiabatic). The 
corresponding blast wave Lorentz factor is shown in Fig. [6] and the 
expected bolometric light curve is presented in Fig. [7] In this case, 
the luminosity after the peak initially decreases as a steep function 
of time, and then it follows a flatter decay. 

• Adiabatic losses — Both electrons and protons suffer from adi- 
abatic losses due to the fireball expansion. Adiabatic losses describe 
the re-conversion of the internal random energy generated at the 
shock into bulk kinetic energy. A proper treatment of the adiabatic 
losses is important in two ways. First, when the velocities are non- 
relativistic we can recover the correct scaling /3 oc R~ 3/2 between the 
fluid velocity and the shock radius predicted by the Sedov-Taylor 
self-similar solution. Models that do not account for the adiabatic 
losses (P99) predict instead a steeper scaling, as /? oc R 3 (see the 
dotted line in the grey area of Fig. [TJ. Second, a proper descrip- 
tion of the adiabatic losses allows to capture the blast-wave re- 
acceleration in the case of an ambient medium where the density 
decreases faster than R 3 , which cannot be properly described nei- 
ther by P99 nor by H99 (Fig.©. 

The model we present is particularly suitable to study the ini- 
tial expansion of GRB external shocks, when a peak (or a change of 
slope) in the light curve is expected to arise. We propose an equa- 
tion to derive the initial Lorentz factor To from the peak time of 
the light curve. Our formula is valid for any value of e and j (see 
Eq.[30]and the normalization factors listed in Tab. [TJ. In the adia- 
batic (e < 1(T 2 ) and homogeneous (s = 0) case we compare our 
formula in Eq. [28] with the equations proposed in the existing liter- 
ature (e.g., M07 and SP99, which differ from each other by a fac- 
tor of 2). We predict a value for To which is intermediate between 
M07 and SP99, yet showing a better agreement with the formula of 
SP99. The equation by M07, often adopted in the literature, tends 
to overestimate To by a factor of ~ 1.6. As an example, for the 
two GRBs studied in M07, we find r ~ 250 instead of r ~ 400. 
The formula derived in this paper is consistent with that used in 
Ghirlanda et al. 2012 for the estimate of To for GRBs with a peak 
in their optical light curve. 

However, some caveats should be pointed out. Our estimate 
of To, similarly to the formulae by SP99 and M07, is based on the 
bolometric light curve, while the observations are performed in a 
given frequency band. Moreover, our estimate for the photon ar- 
rival time (see Eq. [26] which is used to compute the luminosity) 
relies on the simplifying assumption of emission from a spherical 
shell. Effects related to the shape of the radiation spectrum, the ge- 
ometry of the emitting material and the viewing angle (in the case 
of a collimated outflow) can affect the relation between the peak 
time of the light curve observed at a given frequency and r . How- 
ever, we point out that the approach by M07 (commonly adopted 
in the literature) tends to overestimate To, independently from the 
caveats listed above, since the difference mostly arises from their 
overestimate of the deceleration radius, as shown in Fig. [8] 

A complete modeling of the afterglow flux in a given fre- 
quency band cannot be presented here, since it requires a descrip- 
tion of the radiative processes responsible for the observed emis- 
sion, which we have not considered in this work. The dynamical 
model presented here coupled to a radiative model accounting for 
synchrotron and inverse Compton losses will allow us to compute 
spectra and light curves at different frequencies and to investigate 



several issues related to the afterglow physics. Moreover, the in- 
troduction of a radiative model will allow to overcome some of the 
simplifications adopted in this work, such as the requirement of fast 
cooling, implicitly assumed in deriving the formula for the lumi- 
nosity in §4. 1 . In particular, we suggest that our dynamical model 
will be particularly suitable to study: 

• the relation between GeV emission and external shocks. In a 
forthcoming paper, we plan to investigate whether the high-energy 
GeV emission in GRBs can originate from afterglow radiation in 
external shocks. In particular, thanks to the formalism developed in 
this work, we will be able to account for the presence of electron- 
positron pairs, the role of radiative losses, and the temporal evolu- 
tion of the radiative efficiency from a highly radiative phase to the 
standard adiabatic regime. 

• complex circum-burst density profiles and their effect on the 
light curve, in order to explain the rebrightenings, plateaus and 
bumps sometimes detected in the optical and X-ray data. 



APPENDIX A: RE-ACCELERATION INTO VACUUM 

We show that the Lorentz factor of an adiabatic blast wave scales 
as r oc R 3/4 as the shock reaccelerates into the vacuum. We assume 
that the circum burst mass density has some profile p(R) up to a 
radius R vac (where the blast wave is still ultra-relativistic), and after 
that the density drops to zero. For R > R vac , Eq.[7]reduces to 



dT 
dR 



e $ dR 



(M +m vm )c 2 + E' int '^f 



(Al) 



where m vac = 4^J Q p{r)irdr is the swept-up mass up to ra- 
dius R vac . For an ultra-relativistic adiabatic blast wave in the de- 
celeration phase, the comoving internal energy satisfies E' jm » 
(Mo + m)c 2 , and the same will hold at the early stages of the 
subsequent re-acceleration phase. In addition, we assume that for 
R S; Rvac, most of the post-shock plasma is still highly relativistic 
(i.e., y' ad = p' ad » 1), despite the adiabatic losses. From Eqs. 1 121 
and|13l we obtain 



dE 



dR 



, 1 1 dlogT 

-E'JR) — 

"" \R 3 dR 



(A2) 



which relies on the assumption that the comoving shell thickness 
is ~ R/T also in the re- acceleration ph ase, so that the comoving 
volume is V oc R 3 /T (see lShapirolll980|) . In the limit E' jnt » (M + 
m vac)c 2 , Eq. lAll then reduces to 
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whose solution is T oc if 3 ' 4 . 

One can wonder why we obtain T oc R 3/ A , while the init ial ac- 
celeration of the fireball follows T oc R (e.g.. lMeszarosll2006l) . The 
difference is due to different choices for the comoving volume V. 
During the initial acceleration phase V oc R 3 (i.e., a filled sphere), 
while in the case we are considering here (i.e., forward shocks giv- 
ing rise to the afterglow) we assume V oc R 3 /T, which corresponds 
to a shell of thickness ~ R/T 2 in the progenitor frame. If we were 
to assume V oc R 3 , we would recover T oc R. 
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